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We present a strategy for estimating tiie error of truncated functional flow equations. While 
the basic functional renormalization group equation is exact, approximated solutions by means of 
truncations do not only depend on the choice of the retained information, but also on the precise 
definition of the truncation. Therefore, results depend on specification parameters that can be used 
to quantify the error of a given truncation. We demonstrate this for the BCS-BEC crossover in 
ultracold atoms. Within a simple truncation the precise definition of the frequency dependence of 
the truncated propagator affects the results, indicating a shortcoming of the choice of a frequency 
independent cutoff function. 
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I. INTRODUCTION 

Understanding the behavior of strongly correlated 
quantum many-body systems is an important and chal- 
lenging task in numerous areas of modern physics, rang- 
ing from heavy ion collisions and neutron stars to high 
temperature superconductors and ultracold atomic gases. 
Due to the absence of a small expansion parameter, 
non-perturbative methods have to be employed for the 
description of these systems. The necessity to go be- 
yond perturbative approaches arises generically in many 
physical situations, including quantum gravity, systems 
with largely different length scales, or non-perturbatively 
renormalizable theories. 

Typically, providing error estimates for non- 
perturbative techniques is hard. For instance, the 
accuracy of numerical Monte Carlo simulations may 
be measured by the relative change of observables 
due to enlarging the number of grid points. Similarly, 
including higher order terms in truncations of Func- 
tional Renormalization Group or Schwinger-Dyson 
equations partially provides a notion of convergence 
of the corresponding results. However, such improve- 
ments are numerically costly and it is desirable to find 
methods which allow for error estimates within a given 
approximation scheme. 

A promising direction to quantify experimentally the 
reliability of non-perturbative methods has emerged with 
the rapidly evolving field of ultracold quantum gases. 
The latter can be realized and controlled in experiment 
with unprecedented accuracy. Observables such as bind- 
ing energies, the equation of state, the phase diagram or 
collective modes can be measured over a wide range of 
interaction strengths, temperatures and densities |TH3]. 
The precise knowledge of the microscopic Hamiltonian 
of ultracold alkali gases, combined with the universality 
of long range physics, allows for a comparison with first 
principle methods for computing the partition function 
which results from this Hamiltonian, and which deter- 
mines the macroscopic physics of the system. This possi- 
bility for precise measurements makes a reliable estimate 



of errors for non-perturbative methods even more urgent. 

In this work, we address the estimate of errors for the 
functional renormalization group for the effective average 
action. The simple one loop form of the exact functional 
flow equation [4], together with physically intuitive non- 
perturbative truncation schemes, makes the method par- 
ticularly suitable for a study of non-perturbative prob- 
lems |5Hl2j . In order to have a practical example we dis- 
cuss possibilities for an error estimate for the use of func- 
tional renormalization for the BCS-BEC crossover of non- 
relativistic fermions in three dimensions. The system is 
archetypical for the transformation of physical laws when 
going from the simplicity of the microscopic theory to the 
complexity of the effective many-body system: Quantum 
and thermal fluctuations wash out many details of the mi- 
croscopic laws and new structures such as bound states or 
collective phenomena emerge. In particular, for infinite 
s-wave scattering length between the atoms the so-called 
unitary Fermi gas constitutes a universal strongly cou- 
pled system. Important benchmark observables for the 
latter at zero temperature are the Bertsch parameter, the 
gap parameter and the Tan contact. Their precise deter- 
mination attracts a lot of interest, both experimentally 
and theoretically, and allows for comparison of different 
non-perturbative approaches in the above sense. Our aim 
is to obtain an error estimate for the results of the FRG 
approach to the BCS-BEC crossover for the truncation 
presented in Refs. [T3HT7] . 

This paper is organized as follows. In Sec.|lT]we suggest 
a method to estimate the error of solutions to truncated 
flow equations by means of specification parameters. We 
introduce the microscopic model of ultracold fermions 
in the BCS-BEC crossover and discuss the used trunca- 
tion scheme in Sec. |III[ In Sec. |IV| we define specifica- 
tion parameters by means of finite difference projections 
for the wave function renormalization and gradient co- 
efficient of the inverse boson propagator. In Sec. |V] we 
investigate the behaviour of the equation of state, the 
Bertsch parameter, the gap parameter and the Tan con- 
tact upon variations of the specification parameters. In 
Sec. |VI| we draw conlcusions about the accuracy of the 
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used truncation scheme and discuss how this suggests 
possibe improvements of the latter. App. [A| contains the 
truncation in more detaiL The flow equations for the 
running couphngs are presented in App. [B] The emer- 
gence of universaUty due to fixed point solutions of the 
running couplings and their corresponding initial values 
are discussed in App. [C| 



II. ERROR ESTIMATES FOR FUNCTIONAL 
RENORMALIZATION 

The evolution of the effective average action Ffc with 
the renormalization scale k is governed by an exact flow 
equation [3] 



dkTk = ^STr 



Rk) ^dkRk 



(1) 



Here, rfc[(/3] is a functional of bosonic and/or fermionic 

(2) 

fields (fi, and Tj, [ip] is the second functional derivative 
or inverse propagator in the presence of arbitrary fields. 
The supertrace STr includes a single momentum integral, 
as characteristic for a one loop extremum. The infrared 
cutoff Rk should guarantee that only a small momentum 
range of around contributes to Eq. ([T]), and that 
the r.h.s. is ultraviolet and infrared finite. 

Despite its simple structure, Eq. ([IJ is a complicated 
non-linear functional differential equation. A general 
functional rfc[iy9] involves infinitely many parameters or 
couplings. Truncations reduce this infinite set to a finite 
set of couplings or functions that we may call the flowing 
data g. Eq. ([T]) then translates to a set of flow equations 
given by 



kdkfj = C(5)- 



(2) 



Here, C,{g) are the beta functions or flow generators for 
the flowing data g. Typically, g contains the effective po- 
tential, the propagators and several vertices of the theory. 
To obtain a numerically tractable set of equations, g can 
contain at most a finite number of functions. These can 
be parametrized in different ways, e.g. by a polynomial 
expansion or function values at given arguments. The 
approximation procedure for functional renormalization 
consists in a truncation which specifies the used flowing 
data g, a computation of the corresponding flow genera- 
tors C{g)j and finally a numerical solution of Eq. ([T]) with 
initial data g{A) given at some microscopic scale A. For 

— ^ the infrared cutoff Rk is removed, all fluctuations 
are taken into account, and the data g{k — )■ 0) corre- 
spond to physical n-point functions or similar quantities. 
Benchmark quantities as the Bertsch parameter can be 
extracted from g{k — >■ 0). 

For a given truncated flowing data g, the flow genera- 
tors are not uniquely determined. Indeed, the reduction 
of the exact flow equation to a finite set of equations 
introduces the need for a specification how (^{g) is deter- 
mined in terms of g. This prescription fixes the treatment 



of couplings which are left out in the flowing data. For 
example, the inverse propagator of a classical statisti- 
cal model is a function P{q^) of the squared momentum 
(This holds if the system is invariant under transla- 
tions and rotations.). One may choose to approximate 
the momentum dependence of P{q^) by one simple cou- 
pling, P(q^) — P{0) = Zkq^ ■ The wave function renormal- 
ization Zk depends on k and is part of the flowing data. 
One has to specify, however, how one defines Zk precisely. 
An example is 



P{c^e) - P(0) 



(3) 



(For c — >■ this amounts to Zk — dP / dq^\qi^Q.) The 
flow generator for Zk (given by —rjzk, with rj the anoma- 
lous dimension) depends on the choice of the specification 
parameter c. A more precise form of Eq. ([T]) is then 



kdkg = C{9,c), 



(4) 



where c is a set of specification parameters. 

Besides the specification parameters c, the flow gen- 
erators also depend on further numbers like cutoff pa- 
rameters, bosonization parameters etc. We may collect 
them in a set of flow parameters w which replace c in 
Eq. Q. Physical results cannot depend on the choice 
of parameters w. The cutoff is only a technical device 
and its choice has to drop out in the final results [S]. 
The same holds for the procedure of bosonization |18j . 
for example by the precise implementation (19. ,20] of a 
Hubbard-Stratonovich transformation [3T1 [55] or similar 
concepts for flowing bosonization [53]. The specification 
parameters c are obviously only needed for a given ap- 
proximation scheme and have no meaning for physical 
results. 

In a given truncation, however, the results will depend 
on the choice of w. The variation of results within a 
reasonable range of w therefore provides a simple error 
estimate for the short-coming of a given truncation. In 
addition, stability or fragility of the flow equations with 
respect to changes of w can indicate possible ways to 
improve the truncation scheme. 

Investigations of the dependence of results on cutoff 
parameters have been performed in the past by employ- 
ing a family of cutoff functions [24]. This had led to 
criteria for optimized cutoffs ^ i24n27j . Also the depen- 
dence on free parameters in partial bosonization has been 
studied [H], showing how the "Fierz Ambiguity" in the 
Hubbard-Stratonovich transformation is resolved. The 
influence of the relative cutoff scale on observables in a 
theory of both bosons and fermions has been investigated 
in [28 30^. 

In this paper we concentrate on the influence of the 
specification parameter c on the ERG approach to the 
BCS-BEC crossover in the truncation of Refs. [TSUlSj . 
For this quantum statistical system the inverse boson 
propagator -P,/,((Zoi 9^) is a function of two variables, the 
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squared space-momentum and the time component go 
that is related to the Matsubara frequency. The trunca- 
tion of Ffe uses for each k a simple form for the inverse 
propagator, 



(5) 



with flowing data comprising the two "couplings" Zk and 
Ak- In Refs. [13 15 the corresponding flow equations for 
the running couplings Zk and Ak have been defined by 
the derivative of the flow equation for the boson propa- 
gator at vanishing momentum and vanishing frequency. 
In this work, we generalize this prescription by means of 
a finite difference projection according to 



dkZk = 



dkP^.k{clk^O)-dkP^,kiO,0) 



(6) 



and similarly for Ak with a parameter Ca- Although the 
flowing data g ~ (Z, A) remains the same, the flow gener- 
ators additionally depend on the specification parameters 
Cz and Ca- Indeed, an exact flow equation for i-0,fe((Zo, q^) 
can be derived from the second functional derivative of 
Eq. ([I]) [5]. Even for a given truncation of the vertices 
appearing in this equation the flow generators for Zk and 
Ak will depend on the point in momentum space where 
dkP(f>.k is evaluated, and therefore on the choice of Cz and 

Ca- 

The motivation for introducing the finite difference 
projection in Eq. ^ resides in the fact that, for a proper 
choice of regulator functions, the loop integral appearing 
in the flow generators at a given scale k is dominated 
by values of the propagator with qo « k^ and q^ ~ k^- 
Hence, the truncation should resolve the inverse prop- 
agator in the vicinity of these in a sufficient manner. 
This is ensured by the finite difference projection with 
Cz and Ca near unity. In many situations a continuation 
to Cz, — induces only a small effect, while values of 
Cz and Ca substantially larger than unity are not mean- 
ingful. 

In the FRG approach to the BCS-BEC crossover, so 
far, only regulator functions which do not depend on 
frequencies have been employed for the sake of techni- 
cal simplicity. Unfortunately, this implies that both very 
large and very small frequencies (as compared to fc^) con- 
tribute to the loop integrals in the flow generators. Only 
the momentum integrals over space-momenta q are effec- 
tively restricted to a small interval. Because of this short- 
coming in the choice of Rk one of the main advantages of 
the flow equation ([T]), namely the effective restriction of 
all momentum integrals to a small range, is lost. Linear 
approximations to the go-dependence of the inverse prop- 
agator P(t,.k{q[),q^) may therefore become insufficient. A 
substantial dependence of results on Cz could reveal such 
an insufficiency, and we indeed will find that this is the 
case. This demonstrates that our approach to an error 
estimate can be used directly to detect particular short- 
comings of a given truncation or choice of cutoff. 



To investigate the relevance of Cz and c^, we compute 
the equation of state, the gap parameter, the Tan contact 
and the dimer-dimer scattering length in the crossover 
as a function of the specification parameters. For sim- 
plicity, we restrict to the case of zero temperature. We 
find that the parameter Ca has no important influence 
on the results within our numerical resolution. In con- 
trast, the results depend significantly on c^. A detailed 
discussion is presented in Sec. [V| and our Conclusions VI 



This result indicates that the regulator works well for 
spatial momenta, but fails to effectively cut off the con- 
tributions of frequencies outside an interval go ~ fc^- As 
a result the detailed frequency dependence of the boson 
propagator matters, and a too simple approximation for 
this frequency dependence leads to substantial inaccu- 
racies of the results. Incorporating a regulator function 
which depends on both spatial momenta q^ and frequen- 
cies go may result in observables which are much more 
insensitive to variations in both Cz and Ca- 



III. MICROSCOPIC MODEL AND 
TRUNCATION 

We consider an ultracold Fermi gas in three dimen- 
sions, where the atoms can occupy two different hyper- 
fine states. We assume a balanced population. The 
atoms are described by a two-component Grassmann field 
-0 — (-01, -02). The collective bosonic degrees of freedom 
are incorporated in a complex scalar field ip- The eu- 
clidean microscopic action is given by 



l-l/T . 
S[4;, ip]^ j dr y d^x 

+ P*{l'K - 2(U)(,C - h{ip*1pl1p2 - 93-0]; 02))- 

(7) 

It is fixed at a sufficiently large momentum scale A, where 
the interactions between the atoms can be approximated 
to be pointlike. Hence, A^^ has to be much larger than 
the characteristic length scale where details of the inter- 
atomic potential are resolved. The latter is of the order 
of the van-der-Waals length or the effective range. 

We work in non-relativistic natural units with h — 
fee = 1 and rescale 2M ~ 1, where M is the mass of 
the fermionic atoms. The fields depend on the euclidean 
time T, which is restricted to a torus with circumference 
1/T, where T is the temperature of the system. We in- 
troduce the Feshbach or Yukawa coupling h, which is 
directly related to the width of the Feshbach resonance 
[311 132j . After renormalization, the term v\ corresponds 
to the detuning from resonance. This parameter can be 
varied experimentally by applying an external magnetic 
field. 

The basic idea underlying the Functional Renormal- 
ization Group is to take the microscopic action S at ul- 
traviolet scale A as a starting point and to successively 
include quantum and thermal fluctuations of momenta 
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larger than a flowing momentum scale k. For fc = 0, 
all fluctuations are included and one arrives at the full 
efi^ective action F. The latter, in turn, is the generat- 
ing functional of the one-particle irreducible correlation 
functions of the theory, and thus contains all information 
about the macroscopic system. At intermediate values of 
k, the effective average action F^, can be interpreted as 
the effective action for a theory valid at momentum scale 
k. In summary, the functional F^ satisfles Fa — S and 
Fq = F. Its evolution is determined by the exact flow 
equation ([T]). More precisely, in this equation the "su- 
pertrace" STr denotes an integration over momenta and 
a summation over fleld indices, with the characteristic 
minus sign for fcrmions. The regulator has to satisfy 



lim i?fe(Q) = fc2, 
hm i?fe(Q)=0, 



(8) 



with Q = {qQ,q). In principle, one would like to have a 
similar property for the dependence of Rk on qq , but this 
is not realized for the cutoff functions depending only on 
that have been employed so far. 

In this work, we use a basic and physically motivated 
truncation for the effective average action of the BCS- 
BEC crossover. The momentum dependence of the in- 
verse boson propagator is assumed to be of the form 



P^MQ) ^ ^ZkQo + AkQ^ /2. 



(9) 



2fc2)2 and 



This ansatz constitutes a systematic infrared expansion 
of the most general form of the bosonic self-energy, see 
Eq. ( A3 ). A Taylor expansion around the point (gq, q^) 
{elk'', C^k"^) would lead to corrections ~ {qo — c 
{q^ — c^k'^Y that we neglect here. This is well motivated 
if the coefficients of these corrections are not too large 
and if the effective momentum range contributing to the 
flow is restricted to the neighborhood of k"^ . In contrast, 
if the effective momentum range is large, as for go in 
case of an gg-independent cutoff, the deviations of the 
true propagator from the ansatz ([9]) may have important 
effects. 

We neglect fluctuation effects on the fermion propaga- 
tor and the Feshbach coupling. The running couplings 
Zk and Ak, respectively, are referred to as wave function 
renormalization and gradient coefficient in what follows. 
We include a scale-dependent effective potential Uu{p, fi), 
which is a function of the f7(l)-invariant p = (f>*(j) and the 
chemical potential /i. It is expanded around its scale- 
dependent minimum po — po.fc and the chemical poten- 
tial. Our truncation is discussed in more detail in App.fX] 

The choice of regulators is a central issue for the so- 
lution of the flow equation. To implement the idea of 
a momentum-shell integration at each renormalization 
group step, one should use cutoff functions Rk{Q) which 
arc localized at q^ k, k"^ and go ~ for a given scale k. In 
practical applications of the FRG, however, one often re- 
lies on regulators which only depend on spatial momenta 



q^ . This simplifles the computation of the flow genera- 
tors, since Matsubara summations over frequencies can 
be carried out analytically. Moreover, this choice typi- 
cally provides qualitatively convincing results. Although 
optimization criteria for regulators can be derived, it is, 
at present, not fully understood how much the results 
of FRG calculations are influenced by the choice of the 
function Rk{Q). 

In this work, we use optimized regulator functions 



R^,k{Q) = Afc(fc2 _ q^l2)e{e - q^/2) (10) 



and 



i?.0,fe(Q) = {khgn{q^ ~p)- (g2 - ^))e{e - |g2 _ ^,|), 

(11) 

which do not depend on go- For the bosons, we sup- 
press fluctuations with small momenta, whereas for the 
fermions we regularize around the Fermi surface. 



IV. FINITE DIFFERENCE PROJECTION 

In addition to choosing an ansatz for the effective av- 
erage action, one has to specify projection prescriptions 
for the scale-dependent parameters. Here, we use estab- 
lished projection prescriptions for the running couplings 
contained in the effective potential Uk{p, A*), see App. [A] 
The wave function renormalization Z], and the gradient 
coefficient Ak are projected by means of flnite differences 
of the inverse boson propagator. 



Zk = - 



Po 



Ak = 2- 



GI,L(0,P')-GI^22(0,0) 



J>=CAk 



(12) 
(13) 



Herein, the inverse propagator matrix Gi ^ is deflned as 



5cj),{P)54>,{P'y 



(14) 



We work in the real fleld basis for the bosons with (pi and 
4>2 given by 



(15) 



In Eqs. (12) and (13), we introduce the speciflcation pa- 



rameters Cz and Ca as a measure of the relative width 
of the corresponding flnite difference with respect to the 
flow parameter k. In the limit Cz,Ca — 0, the projec- 



tions in Eqs. (12) and (13) reduce to derivatives eval- 



uated at vanishing frequency and momentum. This is 
an established projection prescription used frequently in 
the literature, in particular for the FRG approach to the 
BCS-BEC crossover. We refer to it as the "derivative 
projection". While the latter approximates the inverse 
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FIG. 1: Real part of the inverse boson propagator for fixed 
po as a function of momentum for some finite k in arbitrary 
units. The black solid line represents the full inverse boson 
propagator, as computed from an integration over the fermion 
loop. It shows a typical quadratic dependence for small mo- 
menta p/k <^ 1 and flattening for larger momenta. The blue 
dot-dashed and red dashed lines show the approximations for 
the derivative projection and the finite difference projection 
for Ca = 1, respectively. While the former is appropriate for 
p/k <^ 1, the latter gives a more accurate approximation for 
p/k « Ca- 



boson propagator around pq = ~ 0, the finite differ- 
ence takes into account information at finite frequencies 
and momenta. It is thus supposed to give more accurate 
approximations for frequencies and momenta po sa fc^ 
and « fc^, respectively. Fig. [l] iUustrates this idea for 
the real part of the inverse boson propagator. 

Differentiating both sides of Eqs. and ( 13 1, we find 



dkZk — — - 



Po 



dkAk 



5;tG^.22(0,p')-9feG^,22(0,0) 



pi 



p=CAk 



(16) 



(17) 



Note that one also has to take into account derivatives of 
the arguments and the denominator. In App. [B]we show 
that these additional terms cancel within the used trun- 
cation, however. Due to the bosonic regulator function 
used in Eq. (10), only momenta up to y/2k contribute 
to bosonic loops. It thus seems reasonable to choose 

^ V^- We consider values of Cz and Ca of order unity 
in the following. 

From a technical point of view, the task of solving the 
flow equations for the scale-dependent parameters of the 
truncation becomes more involved when introducing the 
finite difference projections defined in Eqs. (16 1 and (17 1. 



While all loop integrals can be carried out analytically 
for the derivative projection, this is not the case for the 
one entering the flow of Ak with the finite difference pro- 
jection. A two-dimensional momentum integral remains, 
which has to be computed numerically at each step dur- 
ing the flow. 



V. RESULTS 

By investigating the dependence of observables on Cz 
and Ca, we are able to estimate the potential error due 
to a limited momentum and frequency resolution of the 
inverse boson propagator. To leading order, the bosons 
are non-dynamical degrees of freedom in the BCS limit. 
Therefore, we expect that changing the projections for 
the dynamical components of the inverse boson propa- 
gator has only a small influence in this regime. For the 
unitary and EEC regimes, there is no a priori expectation 
for the behavior of physical observables. 

First, we consider the equation of state in the form of 
the functional dependence of the chemical potential on 
the density. The black solid line in Fig. [2] shows the re- 
sult of the derivative expansion l l3H17j . i.e. Cz = Ca = 0, 
for r = 0. We subtract half the binding energy of a dimer 
to obtain a positive number. In the BCS regime, we find 
good agreement with the expected limit ^{n) — eF{n), 
where ep = kp = (Btt^ti)^/'^ is the Fermi energy of a non- 
interacting ideal Fermi gas. The derivative expansion 
also reproduces the Lee~Huang-Yang equation of state 
of a weakly interacting gas of bosonic molecules in the 
EEC limit. However, it lacks quantitative accuracy in the 
unitary regime. For the Eertsch parameter ^ = I^/^f 
the unitary point, one obtains — 0.55. (The subscript 
indicates the used derivative projection.) This value ex- 
ceeds the experimental finding ^oxp — 0.370(5) (8) 3, 33J. 

The finite difference projection allows us to investigate 
the sensitivity of the equation of state with respect to 
variations of Cz and Ca- Fig. [2] also shows the equation 
of state at zero temperature for different values of Cz 
and Ca- The plot agrees with our expectation that the 
system does not depend on Cz and Ca in the ECS regime. 
Moreover, we observe the result not to depend on Ca 
for the whole crossover, both for vanishing and nonzero 
Cz. In contrast, there is a significant decrease of the 
equation of state for increasing Cz in the unitary and 
EEC regimes. The good agreement with Lee-Huang- 
Yang theory is spoiled for > 0. Roughly speaking, for 
increasing Cz , the result for the equation of state remains 
unchanged in the ECS regime, improves in the unitary 
regime, and worsens in the EEC regime. 

Table |l] lists the Eertsch parameter obtained from 
varying Cz and in the sets {0,1,1.4,2.0} and 
{0,1,1.4,2.0,3.0}, respectively. We recall that, due to 
our choice of regulators, spatial momenta are naturally 
restricted to < 2k^. Therefore, Ca = 1-4 provides 
a reasonable upper bound for this specification param- 
eter. We observe a negligible dependence on variations 
of Ca, but a significant decrease for variations of Cz- For 
Cz = 1.4 we find ^(c^ = 1.4) « 0.46 independent of c^. 

We next discuss variations of Cz for fixed Ca — 0. 
This is motivated by the small dependence of observables 
on Ca and the numerical simplicity of choosing Ca = 0. 
We compute the Eertsch parameter ^, the gap param- 
eter A/ep and the contact C/fcp for the unitary Fermi 
gas. Fig. [3] shows the numerical results for these quan- 
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cz~0, c^=0 
cz~0, c^ — l 
cz=0, CA = 1.4 
cz = 1.4, CA^O 
cz = 1.4, CA = 1 
cz = 1.4, CA = 1.4 



Ca\cz 





1 


1.4 


2.0 





0.553 


0.530 


0.462 


0.252 


1 


0.551 


0.527 


0.456 


0.247 


1.4 


0.547 


0.522 


0.451 


0.252 


2.0 


0.544 


0.519 


0.449 


0.259 


3.0 


0.546 


0.522 


0.455 


0.277 



TABLE L The Bertsch parameter ^ — fi/ep for several values 
of Cz and Ca- Within the numerical uncertainty, the Bertsch 
parameter is not affected by variations of ca- However, it 
decreases significantly for increasing Cz- 



(fcpa) 



VI. CONCLUSION 



FIG. 2: The equation of state of the BCS-BEC crossover at 
zero temperature for specification parameters Cz — (upper 
curves) and Cz = 1.4 (lower curves). For each choice of Cz, we 
vary Ca according to Ca £ {0, 1, 1.4}. While the sensitivity of 
the truncation with respect to Ca is neghgible, it significantly 
depends on Cz- The curve with Cz = = corresponds to 
the result obtained with the derivative projection in Refs. [13]- 



titles as a function of Cz- AH three quantities show a 
smaU dependence for < 1. They approach the val- 
ues obtained with the derivative projection in the hmit 
Cz — >■ 0. While ^ and A/ep decrease monotonically for 
Cz ?J 1, the contact C/kp has a maximum for Cz ~ 1.7 
and decreases monotonically for Cz > 1.7. The deriva- 
tive expansion thus constitutes an upper limit for ^ and 
A/ep, whereas the contact only shows slight modifica- 
tions. We have indicated in Fig. |3] the value Cz = 1 
which may be considered as a reasonable value for a 
truncation at a given scale k. Comparing the values 
C = 0.55, A/ep = 0.60 and C/4 0.11 at Cz = 
with the values ^ = 0.53 (0.46), A/ep = 0.59 (0.55) and 
C/k^ = 0.11 (0.11) at 4 = 1 (4 = 2) gives an estimate 
of the error due to the limited frequency resolution of the 
boson propagator for the used truncation and regulariza- 
tion scheme. 

Finally, we compute the ratio of the bosonic dimer- 
dimer scattering length ad and fermionic atom-atom 
scattering length a. The former is obtained in vac- 
uum in the BEC limit and satisfies = Afc=o/47r |15j . 
Fig. |4] shows flD/a as a function of Cz- We find a 
monotonous decrease. While the derivative projection 
results in od/o — 0.73, the exact value of od/o = 0.6 
[33] is matched for Cz ~ 0.90. However, we do not ex- 
pect our truncation to reproduce the correct quantum 
mechanical value. In order to achieve this, one has to 
include more vertices in the effective average action, see 
Ref. (35J. 



In this work, we have presented a way of quantifying 
the error of truncated solutions to exact fiow equations 
by means of specification parameters. The latter appear 
as free parameters in the beta functions of the truncated 
theory. By varying them within a reasonable range, one 
can estimate the reliability of a given approximation to 
the full effective action. A strong dependence of observ- 
ables on such variations indicates a shortcoming of the 
truncation, whereas stability hints on a good incorpo- 
ration of the corresponding physical effect. Both cases 
have been exemplified here for the BCS-BEC crossover 
of ultracold two-component fermions. 

The interaction between fermionic atoms in our ap- 
proach is described by the exchange of a bosonic dimer. 

The propagator of the latter field replaces the effective 
four-fermion coupling of the purely fermionic theory. 
Hence, in order to sufficiently resolve this interaction, a 
good approximation to the frequency- and momentum- 
dependence of the boson propagator is mandatory. In 
a simple truncation, a scale-dependent derivative expan- 
sion Pffj.kiQ) = iZkqo + Akq^ /2, together with a constant 
part m?^ f. evaluated for Q = 0, comprises the leading 
terms of a systematic low-energy expansion of the most 
general form of the boson self-energy. This ansatz for 
the inverse boson propagator, however, still requires a 
specification of how the running couplings Z^. and A^. 
are related to the exact fiow equation for the effective 
average action. 

Starting from the exact ffow equation for the inverse 
boson propagator, a possible way of projecting onto the 
fiow of Zk and Ak consists in taking the derivative with 
respect to iqo and q^, respectively, at vanishing exter- 
nal momentum. This procedure has been applied in the 
literature and results in particular predictions for macro- 
scopic observables like the equation of state or the gap 
parameter. To obtain an error estimate for the latter 
quantities, we continuously deform the differentiation by 
means of a finite difference formula of width c to project 
onto the couplings Zk and Ak- 
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- A/ep 

- C 

- C/k^ 













FIG. 3: Bertsch parameter ^ (black solid line), gap A/ep (red 
dashed line) and contact C/fcp (blue dot-dashed line) as a 
function of Cz for fixed Ca = for the unitary Fermi gas. We 
observe a monotonous behavior of ^ and A/ep . The contact 
C/kp has a maximum for Cz ~ 1.7. 



More precisely, we define and by finite dif- 
ferences with different widtlis, Cz and Ca, respectively. 
Distinguishing temporal from spatial degrees of freedom 
is, firstly, motivated from the generic form of the non- 
relativistic dispersion relation. Secondly, due to our par- 
ticular choice of cutoff functions which only regulate spa- 
tial momenta, the frequency dependence of the boson 
propagator is expected to be less accurate. 

From an investigation of the equation of state, the gap 
parameter, the Tan contact and the dimer-dimer scat- 
tering length in the zero temperature crossover, we find 
that the specification parameter Ca has no significant in- 
fluence on the results. The stability of the observables 
upon variations with respect to c^, which has been ver- 
ified here for < < 3, is not clear a priori. Due to 
our choice of the regulator, momentum integrals are lim- 
ited to q'^ < 2fc^, thus Ca is limited by Ca < 1-4. For 
this reason, our upper bound < 3 provides a solid ba- 
sis to conclude that the cutoff works well in the spatial 
domain and that the resolution for spatial momenta is 
presumably sufficient. 

In contrast to the stability upon variations with respect 
to Ca , the specification parameter Cz substantially modi- 
fies the predictions for the above observables. This is an 
indication for the error due to an insufficient frequency 
resolution of the bosonic propagator. For the cutoff used 
here the error in the unitary regime is quite substantial, 
typically 10 — 20%. A better cutoff function, which also 
acts as a cutoff for frequencies, may be expected to re- 
duce the error substantially, even within the given simple 
truncation Thus our example shows how specifica- 
tion parameters can help to obtain an estimate of errors 
within a given truncation and how to find weak points of 
an approximation scheme at hand. 

We have concentrated here on the specification param- 
eters for the boson propagator. Similar specification pa- 



0.8 




FIG. 4: The ratio of the dimer-dimer and fermion scattering 
length Od and a, respectively, as a function of Cz for fixed 
Ca ~ 0. The curve monotonically decreases. 

rameters should be investigated for the fermion propa- 
gator. Also the precise definition of the density matters. 
Investigating a whole family of possible density defini- 
tions will introduce additional specification parameters. 
In principle, this continues to the specification of vertices. 
One expects, however, the main uncertainties to arrise 
from the lowest n-point functions for which truncations 
are used. The dependence of results on a manageable set 
of flowing parameters can then be used for an estimate 
of the error of a given truncation of the flow equations. 
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Appendix A: Truncation 

The fiow equation for the effective average action con- 
stitutes an exact relation. For its practical solution, how- 
ever, we have to employ a truncation of the most gen- 
eral form of the effective average action, which reduces 
the infinite number of fiowing correlation functions to a 
finite, computationally manageable set. Several physi- 
cally motivated choices are possible. Here we use a scale- 
dependent derivative expansion, where, for each k sepa- 
rately, the functional F^ is expanded in powers of gra- 
dients of the field variables. Typically, the long-range 
physics are well-captured by this approximation. 

Due to fiuctuations, the pointlike (thus momentum- 
independent) interactions of the microscopic theory are 
replaced by effective vertices with a non-trivial momen- 
tum dependence. For the case of the BGS-BEG crossover, 
an effective four-fermion interaction A^^fc((5i, (52j Qs); a 
momentum-resolved Feshbach coupling hk{Qi,Q2), and 
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an infinite tower of dimer-dimer and dimer-atom vertices 
are generated during the flow. In principle, the Func- 
tional Renormalization Group allows for incorporating 
these effects. For the low energy physics, a reasonable 
choice consists in approximating this complicated vertex 
structure by pointlike interactions for suitable collective 
degrees of freedom. 

In this work, we model the atom-atom interactions by 
an effective boson exchange. The boson can be associ- 
ated to a dimer or correlated atom pair. Higher order 
scattering between atoms is partially incorporated in the 
effective potential, which describes interactions between 
the dimers. We neglect the regeneration of a four-fermion 
vertex beyond the dimer exchange and approximate the 
Feshbach coupling to be given by its microscopic value 
/lA. For more elaborate truncations of the effective av- 
erage action for the BCS-BEC crossover, including the 
four-fermion and atom-dimer vertices, see Refs. [351 136j. 

Our ansatz for the effective average action expressed 
in terms of the mean fields Cp = ((^) and V' = (V') reads 

r[(^, V^] = ^ (V^^-P^,fc(5r, -iV)V; + rhAdr. -iV)(^ 

■(AI) 

The effective average potential is a function of the C/(l)- 
invariant p = (p*(p. For a given scale k, the inverse 
fermion and boson propagators, respectively, are given by 
their classical parts contained in the microscopic action, 
and fc-dependent self-energy corrections I]fc((5), which 
are functions of both frequencies and momenta. Thus 
we have 

P,p,k{Q) = + 'f' - + %AQ)^ ( A2) 

P^,fc(Q) = %,k{Q). (A3) 

In addition to P^^^ the inverse propagator for the bo- 
son contains a momentum- and frequency-independent 
part from the second derivative of the effective potential 
Uk- For the present work, we restrict to a leading order 
derivative expansion of the boson self-energy, and neglect 
corrections to the classical fermion propagator. We write 



Pi,.k{Q) = 190 



P4>AQ) = '^ZkQa + 2^kq^- 



(A4) 
(A5) 



We emphasize that the fc-dependent derivative expan- 



sion of the boson propagator in Eq. ( A5 1 goes beyond 
the classical quadratic dispersion relation. Indeed, for 
a sufficiently local regulator R^^k{Q), the contributions 
to the fiow of Zk and Ak are dominated by fluctuations 
with qq « fc^ and qK. k. For k <C IQI the flow of P^.kiQ) 
is suppressed by inverse powers of IQj such that Zk and 



Ak in Eq. (A5) can be replaced by Z{Q) 



Ak=\Q\ in a rough approximation. 



k=\Q\ 



and 



A{Q) 

Therefore, the ansatz iZkQo + Akq^ /2 can model very 
general structures of the type P{Q) = ^f{qQ)+g{q'^), with 



arbitrary functions / and g. A disadvantage of our choice 
of regulators is that we do not regularize frequencies and, 
thus, all frequencies qo contribute to the flow at each scale 
fc. It is one of the purposes of this work to estimate the 
resulting error of this regularization scheme. 

We can rescale the inverse boson propagator in 



Eq. (|A5|) by Aj. ^ , such that the gradient term is given 



by q'^/2 for each k. This is achieved by expressing the 
effective average action in terms of the renormalized held 



A, ip. The corresponding renormalized prop- 



agator reads 



Ak 



(A6) 



with Sk = Zk/Ak- Due to this modiflcation, the gradi- 
ent coefficient Ak drops out of the fiow equations. The 
flow of Ak enters the flow of renormalized vertices and 
propagators only through the anomalous dimension 



1 

Ak 



kdkAk- 



(A7) 



Whereas Zk and Ak show a strong running with fc, the 
quantities Sk and rjk are of order unity. 

So far, we did not specify the form of the effective 
potential Ukip, p)- In principle, its form may be resolved 
by solving a partial differential equation in the variables fc 
and p, which is obtained from the flow equation of for 
constant boson mean field. For the second order phase 
transitions expected here, it is qualitatively sufficient to 
restrict to a (^"^-approximation of U . More explicitly, we 
write in terms of the renormalized field p — (jfcf): 



A 



Uk{p) = Uk{p) = m^fc(p-po,fc) + ^(p-P0,fc) 
- UkSp. + ak{p - /5o,fe)<^M- 



(A8) 



The minimum oi Uk is a fc-dependent quantity. This 
allows to distinguish symmetric (ordered) and disordered 
regimes of the fiow via po,fc = and po,fc > 0, respectively. 
We have > in the former and m| = in the 

latter case. At zero temperature, the system always ends 
up in the superfluid phase for fc = 0, hence /9o,fc=o > 
and m\ ^.^g = 0. The superfluid gap in our truncation is 



given by 



A, 



Po.k, 



(A9) 



with renormalized Feshbach coupling — h\/Ak- 

The effective potential depends on the chemical poten- 
tial /i. It is related to the pressure P according to 



P{p) = -Uk{po±,p.) 



k=0 



(AIO) 



The density n — dP/dp can be computed by solving the 
flow equation for Uk for two inflnitcsimally close values of 
the chemical potential and a subsequent finite difference 
according to n = AP/A/i. Here we work with a flowing 
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density which satisfies nk=Q = n. Moreover, we approx- 
imate d^Ul. ~ d^rn?^ f, = ak- The flow equations for rife 
and ak are obtained from the flow of the effective poten- 
tial by virtue of a derivative with respect to the offset Sfi. 
The latter measures a difference from the actual chemical 
potential /i which also enters the fermion cutoff. In the 
same spirit we may write dp 



P- Pa,k in Eq. (A8) 



The initial conditions for the running couplings r/k, Sk, 



hi -2 



be chosen such that the boundary condition Fa = 5 of 
the flow equations is satisfied. We show below that only 
the initial condition for the boson detuning z^a constitutes 
a relevant parameter of the system. All other quantities 
are attracted to fixed point values in the early stages of 
the flow, i.e. for large k. For this reason, their initial 
values are not important and we may start directly at 
the fixed point. 



Appendix B: Flow equations 



with respect to p and Sp, we obtain the flow equations 
for m| J,, A^,fe, Uk and ak- 

For the derivation of the flow equations for the wave 
function renormalization Zk and gradient coefficient Ak 
we first consider the inverse boson propagator in the 
{(j>i, (^2)-basis. Within our truncation it is given by 



^,fe, A<^,fe, po,k, Uk, and ak introduced here have to G^^{P) = 



Ak 



ZkPo 



2pU' 



-ZkPQ 



Afe(C 



U'k 



(B6) 



We observe that the fiow equations for Zk and Ak can 
be obtained from the 12- and 22-components of the flow 
equation for G'^^, respectively. (Beyond our truncation 
the 11-component contains an additional term ^ A'i.{pq) 
arising from a possible field dependence of Ak (p) , which 
is not present in the 22-component. Therefore, the latter 
is more suited for the projection of ^fe.) We have 



The flow equations for the running couplings deflned 
in App.|X]are obtained from the flow equation for Ffe by 
appropriate projection prescriptions. The resulting set of 
coupled ordinary differential equations in k can be solved 
by means of standard numerical techniques. 

The flow equation for the effective potential Uk{p) is 
found for a constant bosonic background field </> — ^/2p 
and vanishing fermion mean field. We have 



kdkUkip) = VkpU'kip) + Ui''\p) + Ul'\p). (Bl) 

Herein, the first term takes into account that the effec- 
tive average potential has a trivial running due to the 
fact that its argument p = Akp \s scale-dependent. A 
prime denotes a derivative with respect to the variable p. 
The second and third term, respectively, constitute the 
contribution from bosonic and fermionic fluctuations to 
the flow. Within our truncation, these contributions to 
the flow generator for U{p) read 



Ui'^'ip) 



uPip) 



W3 „ .~JP 
h{P)T2 



(B3) 



with = e(x + l){x + 1)3/2 _ _^ (-l)«)6l(a;)a;3/2 + 

( — l)'^/2 ^j^fj Heaviside step function 9{x). 
We introduced the generalized dimensionless "masses" 



Wi = 



W3 



U'kip) 
hjp 

fc4 ' 



W2 



Ul{p) + 2pU'^{p) 

/fc2 



P 



P_ 
k^' 



(B4) 
(B5) 



These ratios control the decoupling of modes as k is 
evolved towards zero. Taking derivatives of Eq. (Bl) 



dkG^},{P)6iP + P') 



S^dkTk 



,(P')<50,(P) 



(B7) 



We define the finite difference projection 
(po,0)-GTi (0,0) 



Po 



Ak 



,G7i2(0,p2)_G-i (0,0) 



= CAk 



(B8) 



(B9) 



For Cz,Ca — i' this becomes a differentiation at P = 0. 



To relate the finite differences to the flow equation (B7|, 
we compute 



dkAk - 2dk 



G^,22(0.4fe^)-G^.22(0.0) 
ci fc2 



dkG^U0,clk^)-dkG-U0,0) 



c?.k'^ 



(BIO) 



9p2 



i/n 2.2N '^0,22(O>ciA;^) - G^_22(0>0) 



The second term in brackets is in general non-zero, but 
vanishes within our truncation. Thus, the flow equation 
of the flnite difference is given by the finite difference 
of the flow equations. The same holds for dkZk with a 
similar derivation. 

The flow equation for Sk — Zk j Ak is evaluated at van- 
ishing external momentum such that the momentum in- 
tegral can be carried out analytically. Since our choice 
of regulators is frequency independent, the Matsubara 
summation can also be performed explicitly. We find 



kdkSk =r]kSk + Slf^ + 



(Bll) 
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with 



si'' 



cli,kPO,k 



37r2fc(l + u;2)3/2 V 5 
-16 + 2(4^2 _ 4)^^ ^ (28 + 4^2)^2/2 + Q^3 



(4- 



4W2 



hi 



37r2fc(l + u;3)3/2 



(B12) 



(4+4(l + «;3))2 

The function £2{x) has been defined below Eq. (| B3| . In 
the limit — ^ we recover the results from Ref . ;13ftl5j . 

The flow equation of for general is given in terms 
of a two-dimensional momentum integral, which has to be 
evaluated numerically at each step of the integration of 
the flow equations. We have 77fc(cyi) 
with 



(B) 

Vk 



00 A| 



i?0 j L^{q+p) - L^{q) 



JQ [ detB((7)detB((f +p) 



■( 



dets((7;) VdetB((7) detB((f+p). 



4^ 

Jq 



[ A / L^{q + p}- L^-{q) (-. 

Jo n detpiq)detF{q+p)\ 



Act B{q) 



-.c^k^ 

(B13) 



1 



1 



detF(g) VdetF(g) detir(g j ^_ ' 

(B14) 

for the bosonic and fermionic contributions to the anoma- 
lous dimension, respectively. We introduced R^/^ — 
kdkR^/^iq^) and 

I^<I>{q) = + ^0('f ) + ™0 + PoA^, 



detBiq) 
detF(g) 



(B15) 

(B16) 
(B17) 
(B18) 



The Matsubara summations in Eqs. ( B13 1 and ( B14 1 can 
again be evaluated analytically. Since 

Rd 



A,. 



= (2k^ - r7fc(fc2 _ (j2/2))(?(fc2 _ q2/2) (B19) 

depends linearly on 77^, we have to solve the linear sys- 
tem of Eqs. (B13)-(B14) to obtain a closed expression for 



ilk{cA)- The formula for ri}~{cA = 0) is given in Refs. [T5]| . 

The flow equation for the Tan contact has been derived 
in Ref. [23 . At zero temperature it is given by 



kdkCk 



hi 



iVkPo + kdkPo) 



hlk^w. 



24V27r25fc(l + «;2)3/2 



l-^^V (B20) 
5 / 



Appendix C: Universality 

The running couplings of our truncation have to be 
equipped with appropriate initial conditions which en- 
sure that the flow satisfies Ta = 5*. We show here that, 
for a broad Feshbach resonance (as it is realized with 
^Li atoms), only the bosonic detuning term j^a consti- 
tutes a relevant parameter, whereas all other couplings 
lose memory of their initial values. This universal behav- 
ior reflects the fact that for a broad Feshbach resonance 
the microscopic model is completely characterized by the 
value of the scattering length. 

In the early stages of the renormalization group flow, 
i.e. for < A, we have /i ^ fc2 and po.fc = 0. Therefore, 
we can set W2 = wi, W3 = jl = in the flow equations. 
The chemical potential does not influence the flow for 
large k, because the high momentum fluctuations cannot 
resolve its particular value. Hence, the many-body sys- 
tem evolves similar to the vacuum system in this regime. 

The flow equation for the dimensionless renormalized 
Feshbach coupling hi — hi/ k is given by 

kdkhl = - l)hl 



(CI) 



for large k. The anomalous dimension in this regime is 
given by rjk ~ hi with a positive prefactor. Hence, if we 
start with a sufficiently large value for h\ — h\/ K, we 
have dkhl > 0. This reduces the value of hi until we 
have r]k — I. At this point, the flow of the dimensionless 
Feshbach coupling stops and hl/k remains constant un- 
til k becomes of the size of the chemical potential or the 
dimer binding energy (on the BEC side of the crossover) . 
The intermediate regime of the flow, where the anoma- 
lous dimension is given by its fixed point value 



(C2) 



is called the universal regime. In the subsequent para- 
graphs we assume k to be in this range. 



The fine-tuning of the bosonic mass term 



U'j^{0) = v/^ — 2p, has been discussed in detail in Ref. [T5] . 
For every choice of a, there exists a unique initial value 
I'A = f^A(a), which describes a system with fermion scat- 
tering length a. We do not repeat the derivation of this 
result here, because the formulas are not altered by our 
finite difference projection with Cz, 0. 

In contrast to the detuning from resonance, we find a 



modified fixed point structure for the flow of /i^, Sk, A 



and ckfc in the early stages of the flow due to the more 
general projection prescription for Ak and Zk- The flxed 



point of hi is determined by Eq. ( C2 1 with 



■qk 



_hL, 
67r2fc 



(C3) 



The positive constant depends on the choice of the 
specification parameter Ca- We find numerically 



< Xfj{cA) < 1 for all Ca- 



(C4) 
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For = we have a;,, = 1. The scahng behavior of 
the Feshbach couphng in the universal regime is deduced 
from Eq. (ICSl) to be 



Assuming to be at its fixed point value, we find 

-xs = Sk-—, 



kdkSk = thSk 



with fixed point solutions 



hi 
37r2fc 



= Ctk 



2 



(C5) 

(C6) 
(C7) 



(C8) 
(C9) 



Both fixed points correspond to stable solutions. From 
the analytical result for the flow of Sk given in Eq. (B12) 
we read off 



1 



xs 



(1 + 4/4)2 



(CIO) 



We have < < 1, with 2:5 = 1 for Cz — 0. We 
conclude that choosing , 7^ modifies the universal 
values oi Sk and afe. 

The flow of the dimer-dimer interaction strength j, 
for large k receives contributions from both fermionic and 



bosonic diagrams. With m| ^ = fc^^ the corresponding 
flow equation 

ht 72fc5(l-7y,/5) ^2 

(Cll) 



has a flxed point solution X^ ^, ~ X^ -i,(Xfi,xs) found from 
the quadratic equation dkX<ji,k = dkikX^i.k) = 0. 

We have seen that, due to the existence of flxed point 
solutions of the flow, the running couplings rapidly ap- 
proach their universal values. In particular, the universal 
regime is entered far before the chemical potential enters 
the renormalization group flow as a scale. Therefore, ex- 
cept for the mass term, the initial values of the couplings 
at scale k = A are irrelevant. In particular, even for 
initially non-dynamical bosonic degrees of freedom with 
Za = = 0, the corresponding terms in the propaga- 
tor are immediately generated. For this reason, we may 
start with Z\, Aa 7^ as well. To summarize, we choose 



, A = i^A(a) - 2/^, A0,A 



Sa = xs/xj,, UA = 6l(/x)^--^, aA = -2/xrf. (C12) 



3/2 



The initial value of the flowing density results from the 
form of the microscopic fcrmion propagator. 
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